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I.  INTRODUCTION 


The  flight  of  projectiles  covers  a  wide  range  of  speeds.  The  accurate 
prediction  of  projectile  aerodynamics  at  these  speeds  is  of  significant  impor¬ 
tance  in  the  early  design  stage  of  a  projectile.  The  critical  aerodynamic 
behavior  occurs  in  the  transonic  speed  regime,  0.9  <  M  <  1,1  where  the  aerody¬ 
namic  coefficients  have  been  found  to  change  by  as  much  as  100%.  Of  parti¬ 
cular  Interest  is  the  determination  of  the  pitching  moment  coefficient  since 
it  is  used  to  determine  the  static  stability  of  the  projectile.  The  critical 
behavior  in  this  case  is  usually  characterized  by  a  rapid  Increase  in  the 
coefficient  followed  by  a  sharp  drop.  This  rapid  change  in  the  pitching 
moment  coefficient  can  be  attributed  in  part  to  the  complex  flow  structure  and 
in  particular,  to  the  asymmetrically  located  shock  waves  which  exist  on  pro¬ 
jectiles  flying  at  transonic  speeds  at  angle  of  attack.  Computations  of 
three-dimensional  flowfields  at  transonic  speeds  are  thus  needed  to  predict 
the  critical  aerodynamic  behavior. 

In  recent  years  a  considerable  research  effort  has  been  focused  on  the 
development  of  modern  predictive  capabilities  for  determining  projectile  aero¬ 
dynamics.  Numerical  capabilities  have  been  developed  primarily  using  Navier- 
Stokes^"^  computational  techniques  and  used  to  compute  flow  over  slender 
bodies  of  revolution  at  transonic  speeds.  Flowfield  computations  have 
included  both  axi symmetri c**  and  three-dimensional  situations,  initial 

computations^”^  did  not  include  the  wake  or  base  region  of  a  projectile  and, 
thus,  ignored  the  upstream  effect  of  the  base  region  flow  on  the  afterbody. 

An  axi symmetric  base  flow  code**  was  then  developed  to  compute  the  entire 

projectile  flowfield  including  the  base  region.  The  base  flow  code  used  a 
flowfield  segmentation  procedure  which  preserved  the  sharp  base  corner.  This 
technique  was  later  extended®  into  three  dimensions  to  calculate  the  pitch 
plane  aerodynamics  at  transonic  speeds.  Due  to  lack  of  computer  resources, 
only  one  solution  was  obtained  and  reported  in  Reference  5.  In  addition,  the 
calculations  in  References  1,  2,  3  and  5  generally  did  not  have  sufficient 

grid  resolution  due  to  the  lack  of  adequate  computer  resources.  The  avail¬ 

ability  of  supercomputers  such  as  the  Cray  X-MP/48  with  a  128  million  word 
solid  state  disk  device  (SSD)  and  the  Cray  2  allows  for  the  increased  grid 
resolution  which  is  needed  for  accurate  computations  of  three-dimensional 
transoni  c  fl  ows.  ^ 

The  numerical  scheme  plays  an  equally  important  role  for  accurate  predic¬ 
tions  of  transonic  flows.  All  the  calculations  in  Reference  1-5  were  made 
using  the  compressible,  thin-layer  Navi er-Stokes  equations  which  were  solved 
using  the  implicit  Beam  and  Warming  central  finite  difference  scheme. Such 
schemes  require  artificial  dissipation  to  be  added  to  control  numerical  oscil¬ 
lations.  Upwind  schemes  can  have  several  advantages  over  central  difference 
schemes  including  natural  numerical  dissipation  and  better  stability  proper¬ 
ties.  The  numerical  scheme  used  here  is  an  implicit  scheme  based  on  flux- 
splitting^°  and  upwind  spatial  differencing  in  the  streamwise  direction. 

Other  factors  that  have  direct  impact  on  the  3-D  numerical  simulation  are 
the  geometric  complexity  and  efficient  management  of  large  3-D  data  sets. 
These  factors  make  it  necessary  to  develop  zonal  or  patched  methods  where  a 
large  3-D  problem  is  divided  into  a  number  of  smaller  problems.  Each  smaller 
piece  is  then  solved  separately.  The  break-up  of  the  large  data  base  can  be 
achieved  in  various  ways.^^"^**  Reference  11  and  12  are  earlier  applications 
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where  the  data  base  structure  follows  a  pencil  format.  These  numerical  calcu¬ 
lations,  although  promising,  were  based  on  limited  computer  resources.  Refer¬ 
ence  13  shows  the  development  of  a  chimera  grid  scheme.  This  scheme  provides 
multiple  regions  where  communications  between  grids  are  done  by  interpolating 
in  regions  of  overlap.  A  blocked  grid  approach  reported  by  Belk,  et  al^**  does 
not  require  interpolations  at  the  interfaces.  The  schemes  in  References  13 
and  14  are  generally  complicated  since  they  allow  for  embedding  a  block  or 
zone  into  another.  Recently,  a  simple  composite  grid  scheme^^  has  been 
developed  where  a  large  single  grid  was  partitioned  into  smaller  grids.  Each 
of  the  smaller  problems  was  solved  separately  with  simple  data  transfers  at 
the  interfaces.  The  initial  results  obtained  were  very  promising.  The 
present  effort  extends  the  use  of  this  composite  grid  scheme  to  include  the 
correct  modeling  of  the  base  region  of  a  projectile.  Three-dimensional  flow- 
fields  have  been  computed  for  two  different  projectiles  over  a  transonic  speed 
range  of  0.8  <  M  <  1.2.  Computed  results  show  the  critical  aerodynamic 
behavior. 


II.  NUMERICAL  METHOD 


1.  GOVERNING  EQUATIONS 

The  three-dimensional  Navi er-Stokes  conservation  equations  of  mass, 
momentum,  and  energy  can  be  represented  ifi  nondimensional  flux  vector  form  as: 

3,Q  +  +  Fv)  +  +  %)  +  +  Hv)  =  0  (1) 

where  the  independent  variable  t  is  the  time  and  the  spatial  variables  4,  n,  c 
are  chosen  to  map  a  curvilinear  body  conforming  discretization  into  a  uniform 

computational  space.  Here  Q  contains  all  the  dependent  variables  [p,  pu,  pv, 

pw,  e]  and  F,  G  and  H  are  the  inviscid  fluxes.  The  flux  terms  F^,  G^  and  H^ 

contain  the  viscous  derivatives.  The  conservative  form  of  the  equations  is 
used  to  capture  the  Rankine  Hugoniot  shock  jump  relations  as  accurately  as 
possible. 

For  body  conforming  coordinates  and  high  Reynolds  number  flow  where  c  is 
the  coordinate  away  from  the  surface,  the  thin  layer  approximation  can  be  made 
in  the  c  direction.  The  governing  equations  can  then  be  written  as: 

3^Q  +  3^F  +  3^G  +  3j.H  -  Re"^3^S  .  (2) 


Here  the  viscous  terms  in  ;  have  been  collected  into  the  vector  S  and  the 
nondimensional  reciprocal  Reynolds  nixuber  is  extracted  to  indicate  a  viscous 
flux  term.  The  viscous  terms  in  e.  and  n  directions  are  neglected. 

In  differencing  these  equations  it  is  often  advantageous  to  difference 
about  a  known  base  solution  denoted  by  the  subscript  o  as: 


2 


(3) 


^(Q  -  Qo)  ^  «5(F  -  F^)  .  6^(G  -  G^)  .  6^(H  -  H^) 


Re’^6j,(S  -  Sq)  =  -3^Qq  -  3^F.  -  a.G^  -  3,H„  +  Re‘^3,S 


^0  no  i;  0 


5^0 


where  6  indicates  a  general  difference  operator,  and  3  is  the  differential 
operator.  If  the  base  state  can  be  properly  chosen,  the  differenced  quantita- 
tives  can  have  smaller  and  smoother  variation  and  therefore  less  differencing 
error.  The  freestream  is  used  as  a  base  solution  in  the  present  formulation. 

2.  IMPLICIT  FINITE  DIFFERENCE  ALGORITHM 

The  implicit  approximately  factored  scheme  for  the  thin  layer  Navier- 
Stokes  equations  that  uses  central  differencing  in  the  n  and  c  directions  and 
upwinding  in  ^  is  written  in  the  form: 

[I  +  h6^(A'")"  +  h6  6"  -  hRe‘%  -  D.l  ] 

i  c  c  1 '  n 

X  [1  +  h6[(r)"  +  h6^B"  -  D.|^]AQ"  =  -At{6b|(F+)f^  .  ^1;]  +  6^1(r)"  -  F^]  (4) 

+  -  GJ  +  -  HJ  -  Re-%^(S"  -  SJ}  -  D^(q''  -  Q„) 

% 

where  h  =  At  and  the  freestream  base  solution  is  used.  Here  6  is  typically  a 
three  point  second  order  accurate  central  difference  operator,  1  is  the  mid¬ 
point  operator  used  with  the  viscous  terms,  and  the  operators  and  are 

backward  and  forward  three-point  difference  operators.  The  flux  vector  F  has 

been  split  into  F  and  F  ,  according  to  its  positive  and  negative  eigenvalues. 

The  matrices  A,  B,  C  and  M  result  from  local  linearization  of  the  fluxes 
about  the  previous  time  level.  Here  J  denotes  the  Jacobian  of  the  coordinate 
transformation.  Dissipation  operators,  Dg  and  0^  are  used  in  the  central 

space  differencing  directions.®  The  factored  left  hand  side  operators  require 
inversion  of  tridiagonal  matrices  with  5x5  blocks.  This  two  factor  implicit 
scheme  has  been  vectorized  and  can  be  readily  multi-tasked  in  planes  of  ^  = 
constant. 

3.  COMPOSITE  GRID  SCHEME 

In  the  present  work,  a  composite  grid  scheme*®  has  been  used  where  a  large 
single  grid  is  split  into  a  number  of  smaller  grids  so  that  computations  can 
be  performed  on  each  of  these  grids  separately.  Each  of  these  grids  use  the 
available  core  memory  in  turn,  while  the  rest  are  stored  on  an  external  disk 
storage  device  such  as  the  SSD  of  the  Cray  X-MP/48  computer.  The  Cray  2 
supercomputer,  has  the  large  in  core  memory  to  fit  the  large  single  grid. 
However,  for  accurate  geometric  modeling  of  complex  projectile  configurations 
that  include  blunt  noses,  sharp  base  corners  and  base  cavities,  it  is  also 
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desirable  to  split  the  large  data  base  into  a  few  smaller  zones  on  the  Cray  2 
as  wel 1 . 

A  code  developed  for  a  single  grid  can  be  made  to  work  for  a  block  grid 
structure  by:  1)  mapping  and  storing  the  information  for  each  grid  onto  a 
large  memory;  and  2)  supplying  interface  boundary  arrays,  pointers  and 
updating  procedures.  Consider  the  situation  in  Figure  1  in  which  the  single 
grid  from  0  =  1,  is  partitioned  into  four  grids  (or  zones),  G1  through 

G4.  The  base  region  of  the  projectile  is  Included  by  adding  another  zone 
G5.  This  procedure  preserves  the  actual  base  corner.  This  zonal  scheme  has 
been  modified  to  allow  more  than  one  zone  in  the  wake  for  accurate  modeling  of 
other  complicated  base  configurations  including  cavities. 

The  use  of  a  composite  or  blocked  grid  scheme  requires  special  care  in 
storing  and  fetching  the  interface  boundary  data,  i.e.,  the  communication 
between  the  various  zones.  For  the  simple  partitioning  shown  in  Figure  1,  all 
subgrid  points  are  members  of  the  original  grid.  There  is  no  mismatch  of  the 
grid  points  at  the  interface  boundaries  and  no  interpolations  are  required. 
This  procedure  thus,  has  the  advantage  over-patched  or  overset  grid  schemes 
which  do  need  interpolations.  The  partitioned  grid  has  six  interface  bound¬ 
aries,  Ji  =  Ji  ,  Jo  —  1,  Jo  =  Jo  ,  Jo  —  1,  Jo  —  Jo  and  J^  =  1  in  the 
^max  ^  max  max 

streamwise  direction  and  two  interface  boundaries  in  the  normal  direction 
between  grids  64  and  G5.  Data  for  these  planes  are  to  be  supplied  from  the 
other  grids  by  injecting  interior  values  of  the  other  grid  onto  the  interface 
boundaries.  The  details  of  the  data  storage,  transfer  and  other  pertinent 
information  such  as  metric  and  differencing  accuracy  can  be  found  in  Reference 
15,  Further  details  of  differencing  accuracy  are  presented  next. 

The  differencing  accuracy  near  the  interfaces  is  quite  important.  Three 
point  backward  and  forward  difference  operators  are  used  at  the  interior 

points.  Near  the  interface,  for  example,  at  Jp  =  Jo  -  1  a  three  point  for- 

^  ^max 

ward  difference  operator  cannot  be  used  with  one  grid  point  overlap  as  shown 
in  Figure  1.  The  differencing  accuracy  can  be  oropped  from  second  order  to 
first  order;  however,  this  leads  to  inaccuracies  in  the  flowfield  solution 
near  the  interfaces. To  maintain  second  order  accuracy  near  the  interfaces, 

3F 

we  difference,  for  example,  —  at  Jp  =  Jp  -  1  as, 

3^  ‘^max 

—  =  +  3^(F-) 


where  3^  is  the  usual  three  point  backward  difference  operator  and  is  now  a 
**  ^ 

central  difference  operator,  i.e.. 


H 


2a^ 


j-1 


2a^ 
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Near  the  other  Interface  of  grid  2  (J2  -  2),  the  3^  operator  Is  corresponding¬ 
ly  replaced  by  a  central  difference  operator  while  3^  Is  a  usual  three  [lolnt 

forward  difference  operator.  The  planes  Jo  «»  1  and  Jo  »  Jo  are,  of  course, 

^  ^  ^max 

boundaries  for  grid  2  and  get  their  data  from  Interior  flowfleld  solutions 
from  neighboring  grids.  Second  order  accuracy  at  and  near  the  Interfaces  Is 
thus,  maintained.  The  use  of  central  differencing  near  the  Interfaces  has  not 
adversely  affected  the  stability  of  the  scheme. 


III.  MODEL  AND  C0^PUTATI0NAL  GRIDS 

The  first  model  used  for  the  computational  study  presented  here  is  an 
idealization  of  a  realistic  artillery  projectile  geometry.  The  experimental 

model  shown  in  Figure  2  is  a  secant-ogi ve-cyl i nder-boattai 1  (SOCBT)  projec¬ 
tile.  It  consists  of  a  three-caliber  (cne-caliber  =  maximum  body  diameter), 
sharp,  secant-ogive  nose,  a  two-caliber  cylindrical  mid-section  and  a  one- 
caliber  7®  conical  afterbody  or  boattail.  A  similar  model  was  used  for  the 
computational  studies  with  the  only  difference  being  a  five  percent  rounding 
of  the  nose  tip.  The  nose  tip  rounding  was  done  for  computational  efficiency 
and  is  considered  to  have  little  impact  on  the  final  integrated  forces. 
Experimental  pressure  data^^  are  available  for  this  shape  and  were  obtained  in 
the  NASA  Langley  eight  foot  Pressure  Tunnel  using  a  sting  mounted  model.  The 
test  conditions  of  1  atm  supply  pressure  and  320  K  supply  temperature  resulted 
in  a  Reynolds  number  of  4.5  x  10®  based  on  model  length. 

The  physical  grid  used  for  this  computation  is  shown  in  Figure  3.  Figure 
3a  shows  the  longitudinal  cross  section  of  the  3D  grid  while  Figure  3b  shows 
an  expanded  view  of  the  three-dimensional  base  region  grid.  As  shown  in 
Figure  3a,  the  clustering  of  grid  points  near  the  body  surface  is  done  to 
resolve  the  viscous  boundary  layer  near  the  body  surface.  Grid  clustering  has 
also  been  used  in  the  longitudinal  direction  near  the  boattail  and  the  base 
corners  where  large  gradients  in  the  flow  variables  are  expected.  In  addi¬ 
tion,  the  composite  grid  scheme  preserves  the  sharp  base  corner.  The  grid 
consists  of  202  points  in  the  streamwise  direction,  36  points  in  the  circum¬ 
ferential  direction  and  50  points  in  the  normal  direction.  This  amounts  to 
about  16  million  words  of  storage  for  the  code  on  the  Cray  X-MP/48.  Only  up 
to  4  Mw  of  central  core  memory  was  easily  accessible;  therefore,  the  full  grid 
was  partitioned  into  five  smaller  grids  (including  a  base  region  grid)  each  of 
which  would  use  the  core  memory  in  turn  while  the  rest  was  stored  on  the  SSD 
device.  These  computations  were  performed  on  the  Cray  X-MP/48  at  the  US  Army 
Ballistic  Research  Laboratory  (BRL).  Each  numerical  simulation  (includes  all 
partitioned  grids)  took  over  20  hours  of  computer  time, 

A  second  grid,  shown  in  Figure  4,  was  obtained  to  simulate  the  actual 
experimental  SOCBT  configuration  which  includes  a  sting  in  the  base  region. 
This  view  is  again  a  longitudinal  cross-section  of  the  3-D  grid.  The  grid  is 
wrapped  around  the  base  corner  in  this  case.  It  consists  of  238  points  in  the 
axial  direction,  39  points  in  the  circumferential  direction  and  50  points  in 
the  normal  direction.  Computations  on  this  grid  were  performed  on  the  Cray  2 
computer  at  BRL  using  the  same  code.  The  computing  time  for  each  of  these 
simulations  was  20-30  hours. 
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The  second  projectile  under  consideration  is  the  M549  projectile  shown  in 
Figure  5.  This  projectile  has  a  short  boattail  of  about  1/2  a  caliber  in 
length.  For  simplicity,  the  flat  nose  was  again  modeled  with  nose  tip  round¬ 
ing  and  the  rotating  band  was  eliminated.  Aerodynamic  coefficient  data^^'*'^ 
are  available  for  this  configuration  and  were  used  for  comparison.  Computa¬ 
tions  for  this  projectile  have  been  made  for  atmospheric  flight  conditions. 
Figure  6  shows  an  expanded  view  of  the  grid  around  this  projectile  and  shows 
both  the  wind-side  and  lee-side  planes.  The  full  grid  consists  of  298  points 
in  the  axial  direction,  39  points  in  the  circumferential  direction  and  50 
points  in  the  normal  direction.  Calculations  for  this  projectile  were  per¬ 
formed  on  the  Cray  2  computer  at  BRL.  Each  of  these  calculations  took  over  30 
hours  of  computer  time. 


IV.  RESULTS 

The  implicit  time  marching  procedure  was  used  to  obtain  the  desired  steady 
state  result.  Initial  conditions  were  freestream  everywhere  and  the  boundary 
conditions  were  updated  explicitly  at  each  time  step.  The  solution  residual 
dropped  at  least  three  orders  of  magnitude  before  converged  solutions  were 
obtained.  In  addition,  the  surface  pressure  distributions  were  checked  for 
time  invariance.  For  the  computation  of  turbulent  flow,  a  turbulence  model 
must  be  supplied.  In  the  present  calculation,  a  two  layer  algebraic  eddy 
viscosity  model  due  to  Baldwin  and  Lomax^*^  was  used.  Results  are  now  pre¬ 
sented  for  two  cases:  (1)  SOCBT  projectile  with  and  without  sting;  and  (2) 
M549  projecti 1 e. 

1.  SOCBT  PROJECTILE,  0.9  <  <  1.2,  a  =  4 

Results  have  been  obtained  at  various  transonic  speeds  for  both  cases  with 
and  without  modeling  of  the  sting.  Figures  7-10  show  the  Mach  contours  for 
the  projectile  in  the  wind-ward  and  lee-ward  planes.  These  figures  show  the 
expansions  at  the  ogive-cylinder  and  cylinder-boattai  1  comers.  These  figures 
indicate  the  presence  of  shock  waves  on  the  cylinder  and  also  on  the  boattail 
which  typically  occur  on  the  projectile  at  transonic  speeds.  Sharp  shocks  are 
observed  on  the  boattail.  These  boattail  shocks  are  shown  to  be  longitudinal¬ 
ly  asymmetric  due  to  the  influence  of  angle  of  attack.  The  asymmetry  can  also 
be  seen  in  the  wake  flow  behind  the  bluff  base.  As  the  Mach  number  is 
increased  from  0.94  to  0.96  and  then  to  0.98,  the  shocks  become  stronger  and 
move  towards  the  base  of  the  projectile.  At  higher  transonic  speeds  past  the 
Sf^eed  of  sound  (see  Figure  10),  these  shocks  become  weak;  however,  a  bow  shock 
forms  in  front  of  the  nose  of  the  projectile. 

Computations  have  also  been  made  to  investigate  the  effect  of  the  sting  on 
the  transonic  projectile  flowfield.  A  typical  plot  of  Mach  contours  for  this 
simulation  is  shown  in  Figure  11a  for  M^  =  0.96  and  a  =  4°.  As  expected,  the 
sting  has  a  large  effect  on  the  qualitative  features  of  the  flow  field  in  the 
wake  region.  An  experimentally  obtained  shadowgraph  at  the  same  flow  condi¬ 
tions  is  shown  in  Figure  lib.  The  actual  shock  wave  position  is  the  front  of 
the  structure  shown  in  this  shadowgraph.  As  shown  in  Figures  Ha  and  lib,  the 
agreement  between  the  computation  and  experiment  for  the  shock  wave  positions 
is  good.  Figures  12a  and  12b  show  the  wind-side  and  lee-side  velocity  vectors 
in  the  base  region  for  both  configurations.  Figure  12a  is  for  the  crse  with 
no  sting  whereas  Figure  12b  includes  the  sting  in  the  base  region.  In  both 
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cases,  asymmetry  in  the  flowfleld  can  be  observed  between  the  wind-side  and 
lee-side.  Three  pairs  of  separated  flow  bubbles  can  be  seen  in  the  near  wake 
for  the  case  of  no  sting  (Figure  12a).  For  the  case  with  sting  (Figure  12b), 
one  can  see  the  large  primary  bubble  along  with  a  counter  rotating  small 
bubble  near  the  junction  of  the  sting  and  the  base.  The  primary  bubble  is 
more  elongated  on  the  wind-side  and  the  flow  reattaches  further  downstream  of 
the  base. 

Figures  13-15  show  the  surface  pressure  distributions  as  a  function  of  the 
longitudinal  position  and  are  compared  with  experimental  data.^®  Figures  13a 
and  13b  show  the  comparison  at  =  0.96  for  the  wind-side  and  lee-side, 
respectively.  Computed  results  are  shown  for  two  grids,  one  which  wrapped 
around  the  base  corner  and  the  other  which  did  not.  As  shown  in  these 
figures,  the  computed  results  are  virtually  the  same  for  both  computations 

except  near  the  base  corner  where  a  small  difference  can  be  noticed.  The 

agreement  of  computed  surface  pressure  with  experimental  data  is  good  for  the 
wind-side.  On  the  lee-side,  the  agreement  is  good  only  on  the  ogive  nose  and 

small  discrepancy  can  be  seen  on  the  cylinder  as  well  as  on  the  boattail.  The 

expansions  and  recompressions  near  the  ogive-cylinder  and  cyl inder-boattai 1 
junctions  have  been  captured.  Figures  14a  and  14b  show  the  surface  pressure 
distribution  for  =  0.98.  Computed  results  are  shown  for  both  cases  with 
and  without  sting  for  turbulent  flow.  In  the  experiment , the  model  was 
sting  mounted  and  no  boundary  layer  trip  was  used.  Therefore,  it  is  not  clear 
if  the  flow  was  laminar  or  turbulent.  Computed  results  were  obtained  for  the 
sting  mounted  case  for  both  laminar  and  turbulent  flow  conditions  and  compared 
with  experimental  data.  Comparison  of  the  pressure  on  the  wind-side  (Figure 
14a)  shows  generally  good  agreement  of  the  computed  pressure  with  the  experi¬ 
mental  data.  The  largest  differences  between  the  computed  results  are  seen  on 
the  rear  part  of  the  boattail  where  no  experimental  results  are  available. 
The  comparison  on  the  lee-side  (Figure  14b)  again  shows  good  agreement  of  the 
computed  result  with  experimental  data  for  most  of  the  projectile  except  on 
the  second  half  of  the  boattail.  As  expected,  the  computed  result  with  no 
sting  has  the  largest  discrepancy.  Computed  results  with  sting  simulation 
compare  well  with  experimental  data  especially  for  laminar  flow  conditions. 
The  difference  in  computed  pressures  on  the  nose  and  the  cylinder  between  the 
sting  case  and  no  sting  case  is  a  result  of  different  grids  used.  A  typical 
resuit  at  a  high  transonic  speed  =  1.1  is  shown  in  Figure  15.  The  agree¬ 
ment  of  the  computed  surface  pressures  with  experiment  is  very  good.  At  this 
high  transonic  Mach  number,  the  shocks  on  the  cylinder  as  well  as  on  the  boat- 
tail  are  very  weak  as  evidenced  by  the  absence  of  a  sharp  rise  in  pressure  in 
those  areas.  The  expansions  and  recompressions  near  the  ogive-cylinder  and 
cyl i nder-boattai 1  junctions  can  be  clearly  observed  in  Figure  15. 

The  computed  surface  pressures  have  been  integrated  to  obtain  the  aerody¬ 
namic  forces  and  moments.  The  slope  of  the  pitching  moment  coefficient  (C  ) 

a 

is  generally  of  greater  concern  in  projectile  aerodynamics  since  it  is  the 
parameter  that  determines  the  static  stability  of  the  projectile.  Figure  16 
shows  the  variation  of  the  slope  of  the  pitching  moment  coefficient  with  Mach 
number.  It  clearly  shows  the  critical  aerodynamic  behavior  in  the  transonic 
speed  regime,  i.e.,  the  sharp  rise  in  the  coefficient  between  M  =  0.92  and 
0.96  and  its  subsequent  sharp  drop.  This  is  followed  by  a  smooth  decrease  in 
the  coefficienc  as  Mach  number  is  increased  further.  The  increase  in  C 

m 
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between  M  =  0.92  and  0.96  is  of  the  order  of  20%  which  is  a  typical  value 
obtained  from  a  number  of  range  tests  for  similar  projectiles. 

2.  M549  PROJECTILE,  0.7  <  <  1.5,  o  =  2 

Numerical  computations  were  made  for  the  M549  projectile  at  various 
transonic  speeds,  0.7  <  <  1.5  and  at  angle  of  attack,  a  =  2®.  Qualitative 

features  of  the  flowfield  obtained  from  some  of  these  calculations  are  shown 
in  Figures  17-21  where  Mach  number  contours  have  been  plotted  for  M  =  0.85, 
0.90,  0.92,  0.94,  and  0.98  for  both  wind-ward  and  lee-ward  planes.  The  asym¬ 
metry  In  the  wake  region  flow  is  obvious  from  these  figures.  These  figures 
indicate  the  development  and  asymmetric  locations  of  shock  waves  on  the 
projectile  at  transonic  speeds.  At  low  transonic  speeds,  for  e-  ample,  at  M  = 
0.85  (Figure  17)  the  shock  waves  are  just  begnning  to  form  especially  near 
the  boattail  junction.  As  Mach  number  is  inc.  ^-ased  to  M  =  0.90,  the  shocks 
are  already  formed  on  the  projectile  both  near  the  cylinder  as  well  as  boat- 
tail  junctions.  The  flow  expansions  at  these  junctions  can  also  be  clearly 
seen  in  this  figure.  The  small  asymmetry  in  shockwave  locations  can  be 
observed  particularly  with  the  boattail  shocks.  The  wind-side  shockwave  on 
the  boattail  is  a  little  closer  to  the  base  than  its  counterpart  on  the  lee- 
side.  In  addition,  these  shocks  have  moved  downstream  from  the  boattail 
junction.  As  shown  in  Figures  19-21,  with  further  increase  in  Mach  number  to 
0.92,  0,94,  and  0.98,  the  shoe''  waves  (both  on  the  cylinder  and  the  boattail) 
become  stronger  and  gradually  move  downstream.  The  asymmetry  in  the  location 
of  the  shock  waves  become  more  apparent  as  the  Mach  number  is  increased.  It 
*is  this- asymmetry  which  gives  rise  to  the  critical  aerodynamic  behavior.  As 
seen  in  Figure  21  for  M  =  0.98,  the  shock  wave  pattern  is  complicated  and  the 
boattail  shocks  are  located  very  close  to  the  base  corners. 

The  static  aerodynamic  coefficients  have  been  obtained  from  the  computed 
flowfields.  As  pointed  out  earlier,  the  slope  of  moment  coefficient  (C^  )  is 

a 

of  primary  concern.  Figure  22  shows  the  development  of  over  the  projectile 

for  various  transonic  speeds.  Actually,  it  is  the  accumulative  moment  coeffi¬ 
cient  referenced  to  the  nose  and  thus,  the  value  at  the  end  (X/D  =  5.645)  is 
the  final  result.  The  difference  in  this  coefficient  over  the  nose  portion  is 
practically  negligible  for  all  transonic  Mach  numbers.  The  largest  effect  is 
seen  on  the  cylinder  and  boattail  sections.  The  boattail  has  a  dramatic 
effect  as  evidenced  by  the  sharp  rise  in  all  the  curves,  F1gure_23  shows  the 

comparison  between  the  computation  and  the  ava-'lable  data.^^  Here 
"a  a 

is  referenced  to  center  of  gravity  (C.G.)  of  the  projectile.  One  can  clearly 
see  the  sharp  rise  in  between  M  =  0,7  to  0.94  which  Is  followed  by  the 

a 

sharp  drop  with  further  increase  in  Mach  number  in  both  the  computation  and 
the  data.  This  critical  aerodynamic  behavior  observed  in  the  data  is  clearly 
predicted  in  the  computations.  Earlier  data  by  Kline,  et  al^^  is  shown  in 
circles  and  is  a  compilation  of  wind  tunnel  and  free  flight  range  data.  More 
recent  data  for  this  projectile  has  been  obtained  from  BRL  range  firings  by 
Whyte,  et  al^®  and  is  shown  in  squanes.  Each  of  these  data  points  comes  from 
a  six  degree  of  freedom  (DOF)  single  fit  procedure.  The  computed  results 
agree  well  with  the  data  by  Kline,  et  al  for  M  >  1.1.  The  computed  results 
overpredict  Whyte's  data  by  5-7%  at  high  transonic  speeds,  1.0  <  M  <;  1.4. 
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Some  discrepancy  can  also  be  observed  at  low  transonic  speeds  M  <  0,9. 
However,  the  discrepancy  is  again  of  the  order  of  5-7%.  Flight  data  for 
transonic  region  Is  no  better  than  7%^^  due  to  rapid  variation  of  with 

a 

Mach  number.  Thus,  the  computed  results  are  well  within  the  experimental 
accuracy  at  these  transonic  speeds. 


V.  CONCLUDING  REMARKS 

In  conjunction  with  a  new  Navler-Stokes  code,  a  simple  composite  grid 
scheme  has  been  developed.  This  capability  allows  the  use  of  fine  computa¬ 
tional  grids  needed  for  accurate  transonic  flow  computations  on  CRAY  X-MP/48 
or  Cray  2  computers.  The  numerical  method  uses  an  Implicit,  approximately 
factored,  partially  upwind  (flux-split)  algorithm. 

The  three  dimensional  transonic  flowfleld  computations  have  been  made  for 
two  projectiles  for  different  flow  conditions  and  angle  of  attack.  The 
computed  flowflelds  show  the  development  of  the  asymmetrically  located  shock 
waves  on  the  projectile  at  various  transonic  speeds.  For  the  SOCBT  projec¬ 
tile,  computed  surface  pressures  have  been  compared  with  experimental  data  and 
are  found  to  be  1n  good  agreement.  The  slope  of  the  pitching  moment  coeffi¬ 
cient  (C^  ),  determined  from  the  computed  flowflelds,  shows  the  critical  aero- 

dynamic  behavior.  The  computed  for  the  M549  projectile  has  been  compared 

a 

with  available  data.  Again,  the  computed  results  show  the  same  critical 
behavior  seen  in  the  data.  The  discrepancy  between  the  computed  result  and 
data  Is  less  than  seven  percent. 

The  results  of  this  research  provide  the  basis  for  a  new  capability  to 
compute  three  dimensional  transonic  flowflelds  over  projectiles.  This  capa¬ 
bility  In  conjunction  with  the  supercomputers  at  BRL  has  led  to  the  first 
successful  prediction  of  the  critical  aerodynamic  behavior  in  of  artillery 

a 

Shell  at  transonic  speeds.  The  next  step  Is  the  numerical  prediction  of 
Magnus  force  and  moment  for  spinning  projectiles  at  angle  of  attack.  This 
Involves  calculations  of  the  full  three  dimensional  flowfleld  with  no  plane  of 
symmetry. 
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ALL  DIMENSIONS  IN  CALIBERS 


Figure  2.  Model  geometry  of  the  SOCBT  projectile. 
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Figure  5.  M549  projectile. 


Figure  8.  Mach  contours,  SOCBT  projectile,  =  0.96,  a  =  4°. 
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Figure  15. 
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Longitudinal  surface  pressure  distribution, 
SOCBT  projectile,  M_  -  1.1,  a  =  4°. 
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Figure  16.  Slope  of  pitching  moment  coefficient,  C  ,  vs  Mach  number 
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Figure  17.  Mach  contours,  M549  projectile,  M, 
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